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1  INTRODUCTION 


This  is  a  project  to  provide  mission  support  for  the  Communication/Navigation  Outage  Forecast 
System  (C/NOFS)  under  BAA  VS-03-01  during  its  first  four  years  of  operation.  Cornell  will 
support  the  mission  with  ground-based  radar  observations  of  background  ionospheric  parameters 
and  of  equatorial  spread  F  (ESF)  events  from  the  Jicamarca  Radio  Observatory  near  Lima,  Peru. 
Jicamarca  is  capable  of  measuring  ionospheric  electric  field  and  conductivity  profiles  from  the 
valley  region  well  into  the  topside  ionosphere.  Jicamarca  can  also  make  very  detailed  observations 
of  ionospheric  irregularities  associated  with  ESF.  Three  main  tasks  will  be  supported  under  this 
project:  (1)  in-flight  calibration  of  the  electric  field/drift  meter  sensors  on  board  the  spacecraft, 
(2)  provision  of  Jicamarca  radar  support  for  experimental  campaigns,  both  before  the  launch  and 
during  the  satellite  operations,  and  (3)  analysis  and  interpretation  of  the  data  obtained  during 
these  campaigns  in  light  of  the  C/NOFS  mission  goals. 

2  METHODS  AND  PROCEDURES 

The  protracted  pre-launch  period  for  C/NOFS  has  afforded  the  opportunity  to  reexamine  funda¬ 
mental  instability  theory  and  rethink  the  mission  forecast  strategy.  Emphasis  now  lies  in  using  the 
theory  developed  by  [2]  regarding  plasma  shear  instabilities  in  the  postsunset  bottomside  F  region 
ionosphere  as  a  tool  for  understanding  and  forecasting  equatorial  spread  F.  The  basic  idea  is  that 
the  shear  instability  is  more  robust  than  the  ionospheric  interchange  instability  immediately  after 
sunset  and  that  the  former  serves  as  a  seed  for  the  latter,  which  appears  otherwise  to  grow  too 
slowly  to  account  for  the  rapid  production  if  large-scale  irregularities  often  observed.  If  so,  then 
the  parameters  being  measured  and  modeled  by  the  C/NOFS  science  team  should  be  modified 
accordingly.  The  changes  required  are  not  drastic,  as  the  shear  instability  is  closely  related  to  the 
ionospheric  interchange  instability.  Furthermore,  proxies  for  shear  instabilities  should  themselves 
be  monitored.  These  include  so-called  bottom-type  scattering  layers,  which  are  common  features 
of  the  postsunset  equatorial  ionosphere  and  which  are  known  to  become  patchy  very  early  in  the 
evolution  of  large-scale  waves  and  irregularities. 

Bottom-type  layers  represent  a  necessary  condition  for  ESF.  Patchy  bottom-type  layers  more¬ 
over  appear  to  be  a  sufficient  condition.  These  strong  statements  are  based  on  long-term  studies  at 
Jicamarca  and  on  results  from  the  NASA  EQUIS  II  rocket  campaign  on  Kwajalein  in  the  summer  of 
2004.  Patchy  layers  have  been  observed  at  Jicamarca  using  radar  imaging  and  on  Kwajalein  using 
the  scanning  Altair  radar.  Occasionally,  the  patchiness  of  the  layers  is  apparent  in  conventional 
range-time-intensity  (RTI)  radar  data  from  Jicamarca.  Fig.  1  illustrates  a  broken  layer  exhibiting 
~40  min.  periodicity  telltale  of  a  large-scale  precursor  wave  with  a  wavelength  of  the  order  of 
100-200  km.  Fig.  1  shows  another  patchy  bottom-type  layer,  this  time  with  ~  10  min.  periodicity 
indicative  of  a  large-scale  precursor  wave  with  a  wavelength  of  a  few  tens  of  kilometers.  These  two 
scale  sizes  are  prevalent  in  imaging  and  Altair  data. 

Most  often,  however,  RTI  data  from  the  JULIA  radar  at  Jicamarca  (which  has  been  running 
semi-nightly  for  ten  years)  give  no  clear  indication  of  whether  bottom-type  layers  are  patchy  or 
continuous.  Standard  JULIA  modes  were  designed  for  a  synoptic  view  of  E  and  F  region  irreg¬ 
ularities  and  not  explicitly  to  monitor  bottom-type  layers.  We  have  therefore  instituted  a  new 
JULIA  radar  mode  which  will  be  in  operation  during  the  remainder  of  the  C/NOFS  pre-launch 
period.  The  mode  is  optimized  to  detect  irregular  bottom-type  layers  using  two  approaches.  One 
of  these  is  a  split-beam  approach,  whereby  the  echoes  from  two  radar  beams  directed  eastward 
and  westward,  respectively,  will  be  compared.  Significant  differences  in  echo  power  will  indicate 
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Figure  1:  Range  time  intensity  (RTI)  image  depicting  a  patchy  bottom-type  layer  prior  to  the  onset 
of  ESF. 


Figure  2:  Another  RTI  image  showing  a  patchy  bottom-type  layer  in  advance  of  ESF  irregularities. 
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horizontal  inhomogeneity  at  the  scale  size  of  the  beam  separation.  The  other  approach  involves 
interferometry.  Echoes  corresponding  to  one  of  the  beams  will  be  processed  using  spaced  receiver 
correlation.  The  coherence  of  the  signals  will  indicate  horizontal  inhomogeneity  at  the  scale  size 
of  the  radar  beam  width.  An  expert  system  will  assess  the  data  and  eventually  be  applied  to  near 
real-time  forecasting.  Data  collection  using  the  new  mode  is  underway. 

The  best  prospect  for  advance  spread  F  forecasting,  however,  involves  anticipating  ESF  by 
estimating  the  growth  of  the  underlying  shear  instability.  We  can  do  this  by  continuously  evaluating 
the  appropriate  linear  growth  rate  expression,  much  as  the  conventional  strategy  calls  for  evaluating 
the  ionospheric  interchange  instability  growth  rate.  A  combination  of  experimental  and  model  data 
will  be  required,  and  these  must  ultimately  be  available  on  an  operational  basis.  The  existing 
C/NOFS  forecast  strategy  already  incorporates  these  elements. 

We  recently  derived  the  linear  growth  rate  for  the  collisional  shear  instability.  This  is  an 
inherently  nonlocal  instability,  and  the  semi-local  expression  we  found  is  at  best  an  approximation. 
Exact  results  for  a  slab  geometry  have  been  obtained  by  [2].  Our  expression  has  been  validated 
against  numerical  models  and  simulations  [3]. 

This  calculation  takes  places  in  magnetic  dipole  coordinates  (p, <7,<1>)  described  by  [1].  The 
plasma  number  density  and  potential  are  separated  into  zero-  and  first-order  components  according 
to 


4>(p,<7,$)  =  Mp)  +  ut) 

n(p,  q,  <J>)  =  nD(p)  +  ni(p, 


where  <f>0  is  the  potential  associated  with  the  background  zonal  E  x  B  plasma  drift  (v0(p)),  and 
n0(p)  expresses  the  background  plasma  density  gradient. 

To  first  order,  the  continuity  equation  can  be  written  as: 


1  dnQ  k  <f) i 


0 


(1) 


where  h $  is  the  scale  factor  associated  with  the  $  coordinate,  for  example,  and  B  is  the  background 
magnetic  field.  The  dominant  first  order  plasma  drift  velocity  is  the  E  x  B  drift  arising  from  <j>\. 

We  next  invoke  the  quasineutrality  condition,  V  •  J  =  0.  We  evaluate  the  transverse  component 
of  the  divergence  of  the  current  density,  a  quantity  that  vanishes  when  integrated  along  a  magnetic 
field  line  between  points  where  it  exits  the  ionosphere.  The  primary  first  order  currents  we  consider 
are  Pedersen  currents  arising  from  <j)\  and  from  density  perturbations  n\  embedded  in  a  background 
electric  field,  which  in  the  frame  of  reference  of  the  neutral  gas  moving  with  zonal  speed  u  is 
proportional  to  u  —  vQ. 


d_ 

dp 


where  ap(p)  is  the  Pedersen  conductivity.  Combining  (1)  and  (2)  gives: 


(2) 
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=  Vj_  •  J 


The  next  step  in  the  derivation  is  to  multiply  the  preceding  equation  by  <t>\  and  integrate  over 
a  neighborhood  of  p  where  the  perturbed  potential  is  significant.  In  the  formulas  to  follow,  the 
p  integration  is  not  written  but  is  understood.  Applying  integration  by  parts  and  dropping  the 
boundary  terms  then  yields: 


rin  ri(t> 


hphq 


<f> 

u  —  v  o 


K  —  h4>'l  =  Vx  •  j <t>\ 


Hq 

hp  uj  —  {ft/h^Vo  nQ 


where  the  prime  superscript  represent  differentiation  with  respect  to  p. 

At  this  point,  we  make  two  assumptions  which  are  motivated  by  the  exact  solutions  obtained 
computationally  by  [2].  The  first  is  that  the  real  part  of  the  frequency  of  the  growing  wave  is  well 
approximated  by  the  Doppler  shift,  i.e.  u  «  (n/h^Vo  + ry.  The  second  is  that  in  the  vicinity  of 
its  peak,  the  perturbed  potential  varies  approximately  according  to  exp {iKpp).  This  is  essentially 
a  local  approximation  except  that  we  remember  to  maintain  the  p  integration  described  above  in 
the  analysis.  To  avoid  confusion,  we  replace  k  from  above  with  in  equations  that  follow. 

At  this  point,  we  have 


h<j>hq  2 1 X  |2  hphq 


- -irap^2 


hn  U  ~  Vo 

+T~ap~ 

hp 
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*pK*— I4>iI2  =  v±.j# 


(3) 


Integrating  (3)  over  magnetic  field  lines  from  end  to  end  in  the  ionosphere  and  taking  the  potential 
to  be  invariant  along  the  integral  then  produces  the  flux  tube  integrated  growth  rate  estimate  in 
(4),  since  (Vi-J)  =  J\\  I  endpoints  =  0. 


p(u  vo)^o  /no) 

(^ap(Eo/B  +  gJ^n)n0/n0) 

(~hfap> 

Note  that  the  angle  brackets  in  (4)  imply  an  average  over  the  span  of  p  where  the  perturbed  potential 
is  significant  in  addition  to  an  average  over  the  field  line.  This  span  is  likely  to  be  quite  narrow 
in  the  early  stages  of  instability  since  small  spans  are  associated  with  the  largest  growth  rates. 
A  similar  (but  simpler)  approach  yields  the  more  exact  growth  rate  estimate  for  the  generalized 
Rayleigh  Taylor  instability  in  (5),  which  is  a  well-known  result. 

The  growth  rate  estimates  necessitate  the  specification  of  ionospheric  Pedersen  conductivi¬ 
ties,  winds,  and  electric  fields  over  the  magnetic  flux  tubes  of  interest.  Jicamarca  can  measure 
plasma  number  densities  and  electric  fields  overhead,  but  models  are  necessary  to  complete  the 
off-equatorial  specification.  Winds  likewise  require  modeling,  and  we  expect  empirical  models  of 
the  thermospheric  winds  to  improve  once  C/NOFS  wind  measurements  are  available. 

The  estimates  are  no  better  than  the  models  supporting  them,  and  we  have  been  collecting 
baseline  data  at  Jicamarca  in  order  to  help  validate  the  C/NOhS  physics-based  background  model 
(and  also  to  support  the  COSMIC  mission).  Data  from  the  June  campaign  are  shown  in  Fig.  3. 
The  first  day  of  observations  were  marred  by  technical  problems,  and  overall  data  quality  suffers 
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Figure  3-  Jicamarca  1SR  data  from  June.  The  four  rows  depict  plasma  density,  electron  tempera¬ 
ture,  ion  temperature,  and  H+  fraction.  The  He+  fraction  was  not  fit  due  to  the  low  signal-to-noise 

ratio. 
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from  very  low  solar  flux,  ionospheric  content,  and  signal  levels.  Nevertheless,  the  data  highlight  two 
crucial  features  of  the  equatorial  ionosphere  not  well  represented  by  existing  empirical  01  physics- 
based  models.  These  are  1)  a  very  steep  bottomside  density  gradient  after  sunset  with  a  scale 
length  of  the  order  of  10  km  and  2)  a  relatively  dense  postsunset  valley  region,  with  densities  of 
nearly  104  cm  3.  The  one  model  that  seems  to  be  able  to  reproduce  these  features  accurately  is 
the  P1M  model.  This  is  no  accident,  as  P1M  has  been  tuned  to  reproduce  the  climatology  observed 

at  Jicamarca.  . 

More  C/NOFS/COSMIC  validation  campaigns  are  scheduled  for  the  fall  of  2006.  There  is  little 
point  in  attempting  to  forecast  ESF  on  the  basis  of  the  arguments  presented  here  until  contemporary 
physics-based  models  can  be  made  reproduce  data  like  those  in  Fig.  3  and  those  acquired  during 
the  EQU1S  II  campaign  more  accurately.  We  are  working  with  our  colleagues  at  A  FPL  and  on  the 
C/NOFS  EWG  to  improve  the  situation. 
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